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ABSTRACT 



We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb 
with a planct-star mass ratio of q — [9.5 ±2.1] x 10~^ via gravitational mi- 
crolensing. The planetary deviation was detected in real-time thanks to the high 
cadence of the MOA survey, real-time light curve monitoring and intensive follow- 
up observations. A Baycsian analysis returns the stellar mass and distance at 
Ml = 0.64^Q 2g Mq and Di = 5.9^5 4 kpc, respectively, so the mass and separation 
of the planet are Mp — 20^1 and a — 3.3^q ;§ AU, respectively. This discov- 
ery adds another cold Neptune-mass planet to the planetary sample discovered 
by microlensing, which now comprise four cold Neptune/Super-Earths, five gas 
giant planets, and another sub-Saturn mass planet whose nature is unclear. The 
discovery of these ten cold exoplanets by the microlensing method implies that 
the mass ratio function of cold exoplanets scales as dNpi/dlogq oc g-o.7±o.2 ^j^j-^ 
a 95% confidence level upper limit of n < —0.35 (where dNpi/dlogq oc g"). As 
microlensing is most sensitive to planets beyond the snow-line, this implies that 
Neptune-mass planets are at least three times more common than Jupiters in 
this region at the 95% confidence level. 

Subject headings: gravitational lensing, planetary systems 
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Introduction 



S ince the first discov ery of exoplanets orbiting main sequence stars in 1995 (IMayor &: Queloz 



19951 : iMarcy et al.l 120051) . more th an 300 exoplanets have been discovered via the radial ve- 
l ocity method (iMavor et al.l 120041) and m ore than 50 have been detected via their transits 

( Udalski et al. 2004; Konacki et al. 20051). Several planetary can didates have also been de- 

tecte d via direct imaging ( iMarois et al. I boosi : iLagrange et al.ll2009l ) , and astrometry (jPravdo fc ShaklanI 
20091 ). Here, we report the tenth exoplanet discovery by the microlensing method, which is 
another example of a cold, Neptune-mass planet discovered. Although the radial velocity 
and transit discoveries are more numerous, microlensing is uniquely sensitive to these cold 
Neptunes, and the microlensing results to date indicate that this class of planets may be the 
most common type of exoplanet yet discovered. 



LiebesI ( 119641 ) and iMao &: Paczyhskil (jl99l[ ) first proposed exoplanet searches via gravi- 



tational microlensing. T he planet's gravi ty induces small caustics, which can generate small 
deviations in standard ( iPaczyhskil Il986[ ) single lens microlensing light curves. Compared 
to other techniques, m icrolensing is sensitive to smaller planets, down to an Earth mass 
( iBennett fc Rhid 119961 ). and in wider orbits of 1-6 AU. Because microlensing observability 
does not depend on the light from the lens host star, it is sensitive to planets orbiting faint 
host stars like M-dwarfs and even brown dwarfs. Furthermore, it is sensitive to distant host 
stars at several kpc from the Sun, which allows the Galactic distribution of planetary systems 
to be studied. 



discovery (Bond et al. 200^ 


i). So far 8 planetary syst 


ems with 9 planets have been found 


by this technique (1 


Jdalski et al. 2005 


Beaulieu et al. 


200d: 


Gould et al. 2006: Gaudi et al. 


2008: Bennett et al. 


20081: ] 


Done et al. 


2009bl: Janczak et al. 


2010r). which have very distinct 


properties from those detected by other techniques. 


Beaulieu et al. 


(2006) found a ~ 5.5 



Earth-mass planet, which was the lowest-mass planet det ected at that time. This detection 
and the discovery of a slightly more massive planet by iGould et al.l ( l2006l ) demonstrated 
that microlensing is well suited to detecting low-mass planets at orbital distances that are 
currently beyond the reach of other methods. At the time of the discovery of these two cold 
Neptune-mass planets (hereafter "Neptunes") or "Super Earths", two Jovian planets had 
also b een found. These discoveries indicate that cold Neptune in orbits beyond the "snow- 
line" (llda fc Lirul2004j : iLaughlin. Bodenheimer &: Adamsll2004j : Kennedy. Kenyon &: Bromley 
20061 ) around late-type st ars, are significantl y more common than gas giants with frequency of 
>16% at 90% co nfidence (iGould et al.ll2006l ). which is consistent with theoretical simulations 
( Ida fc L"inll2004 ) based on the core accretion model. On the other hand, microlensing has also 



revealed the most massive M-dwarf planetary companion ( lUdalski et al.ll2005t iDong et al. 
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2009ar ) , which would 



ikely be difficuh to form by core accretion ( iLaughhn. Bodenheimer &: Adams 



2004 ). I G audi et ahl (|2008[ ) discovered a system with a Jupiter and a Saturn orbiting an M 



dwarf in a configuration very similar to that of our solar system. Remarkably, this event 
yielded a direct measurement of the masses of the planets and the host star, that was 
confirmed by direct observation of the host star. This system (OGLE-2006-BLG-109Lb,c) 
is the only known multi-planet system with measured masses for the star and planets 
(aside from our own Solar System). The light curve of this event also yielded informa- 
tion about the orbit of th e Saturn- mass planet that confirms that this system is similar to 
ours ( Bennett et al.l 2009). A pla net was also found to orbit a very low-mass host star or 
brown dwarf teennett et ah 2008 ). and this planet was also the lowest mass exoplanet known 
at the time of its discovery. 

Here we report the discovery of another Neptune-mass exoplanet in the microlensing 
event OGLE-2007-BLG-368. We describe the datasets in Section[2l The light curve modehng 
and uncertainty of the parameters are presented in Section [3l and the physical characteriza- 
tion of the lens system is considered in Section [71 In Section IH we discuss the implications 
of microlensing planet discoveries for the exoplanet mass function. The discussion and con- 
clusions are given in Section O 



2. Observations 



The Optical Gravitational Lensing E xperiment (OGLE) (lUdalski 



ing Observations in Astrophysics (MOA) ( iBond et al.ll200ll : ISumi et al. 



2003) and Microlens- 



20031 ) are conducting 



microlensing surveys toward the Galactic bulge to find exoplanets. From 2002 to 2008, the 
OGLE-III survey discovered about 600 microlensing events every year by using 1.3m Warsaw 
telescope with a 0.34 deg^ fie ld-of-view (FOV) mosaic CCD camera with its Early Warning 
System (EWS, lUdalskil |2003| ) . The data have been analyzed in real time and all kind of 
deviations from the usual single lens light cur ves, including planetary anomalies, have been 
detected by the EEWS system Judalskilboosh . The second phase of MOA, MOA-II, carries 
out a very high cadence photometric survey of the Galactic bulge with the 1.8m MOA-II 
telescope with a 2.2 deg^ FOV CCD camera. In 2007, 4.5 deg^ of the central Galactic bulge 
were observed every 10 min, and additional 45 deg^ were observed with a 50 min cadence. 
This strategy enables the detection in real-time of planetary deviations in any of the ~ 500 
microlensing events seen by MOA every year. (Starting in 2010, the new 1.4 deg^ OGLE- 
IV camera will enable OGLE to follow a similar strategy of high cadence monitoring for 
planetary signals.) 



The microlensing event OGLE-2007-BLG-368 was discovered at (R.A., Dec.) (2000) = 



- 7- 



(17:56:25 96, - 32:14:14.7) [(/, 6) = (358.3°, -3.7°)] and alerted by the OGLE EWS system 



flUdalskill2003l ) on 2007 July 10, and independently detected by MOA and alerted as MOA- 



2007-BLG-308 on 2007 July 12. 

Around UT 12:00 20 July (JD= 2454302), MOA observed a series of 9 points that are all 
below the point lens lightcurve, and these are confirmed by a single OGLE point, with higher 
precision. See Figure [H The prompt informa tive data release to t he scientific community 



allowed the SIGNALMEN anomaly detector (IDominik et al.l 120071 ) (now an integral part 



of the Automated R obotic Terrestrial Exoplanet Microlensing Search (ARTEMiS) system 



Dominik et al.l l2008f ) to detect a light curve anomaly that was passed on to 1-3 members 
of each of the major microlensing collaborations, such as the PLANET, /iFUN, RoboNet, 
OGLE and MOA at UT 19:32 20 July (JD=2454302.314), that this was a possible planetary 
anomaly. Given the intensity of microlensing decision-making and the incompleteness of 
the information flow, this distribution proved only partially adequate and failed to reach 
the MOA internal alert system. Based on this alert, the /iFUN SMARTS (CTIO) telescope 
began obtaining data just 5 hours later, shortly after dusk in Chile, after which the PLANET 
Danish (La Silla) and Canopus (Tasmania) telescopes also began observations. Although 
MOA observer did not receive this alert, its high-cadence survey enabled good coverage of 
a steep rise due to the caustic entrance in the light curve the next night, which triggered 
MOA's real-time anomaly alert system to circulate an alert, calling for the firm detection 
of the anomaly, based on its own data at UT 15:58 21 July (JD=2454303. 16528). Here 
the real-time anomaly alert system adds new data points on the lightcurves within 5 min 
after exposures to search for deviations from the single lens lightcurve. This continuous 
early coverage proved crucial for the interpretation of the planetary anomaly. See Figure [H 
Prompted by these anomaly alerts, MOA-II, OGLE-III and other telescopes from PLANET 
and /iFUN began intensive follow-up observations, which densely covered the second peak, 
due to the caustic exit, and less densely for about 50 days. 

Twelve light curves were collected by 7 telescopes. MOA-II 1.8m (Mt. John, New 
Zealand) obtained 1577 images in the MOA-Red wide band, which corresponds roughly to 
a combined I + R filter. OGLE-III 1.3m (Las Campanas, Chile) obtained 12 images in the 
V band and 733 in I. /iFUN SMARTS 1.3m (CTIO, Chile) obtained 22 images in V, 137 
in / and 128 in H. PLANET SAAO Im (SAAO, South Africa) obtained 9 images in V and 
60 in /. PLANET Canopus Im (Tasmania, Australia) obtained 50 images in I. PLANET 
Danish 1.54m (La Silla, Chile) obtained 20 images in V and 129 in /. PLANET OPD/LNA 
0.6m (Brazil) obtained 122 unfiltered images. 

The photometry of this event was much more difficult than the photometry of most 
microlensing events due to a much brighter star located approximately 1.1" to the NW 
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of the source star. This causes very severe problems with stand ard PSF-fitting photometry 
approaches, such as DOPHOT (ISchechter. Mateo, fc Sahalll993f). so the on l y viable approach 



is the Difference Im age Analysis (DIA) method (ITomany fc Crottd Il996l : lAlard &: Lupton 



19981 : lAlard I l2000f ) . The images were reduced by three different implementations of DIA 
photometry. OGLE V and I and CTIO I images were reduced by the standard OGLE 
piA pipeline (lUdalskil 120031 ). Other images wer e reduced by both th e MOA DIA pipeline 



Bond et al.l 120011) and a version of pySIS (v3.0) (lAlbrow et al. 1120091). partly based on ISIS 



( lAlard fc Lupton Ill998l ). but using a numerical kernel (jBramichll2008l ). In the MOA DIA 



pipeline, point spread function (PSF) photometry was performed on the difference images 
with various reference images and PSF fitting radii. 

The resulting MOA DIA light curves, the pySIS light curves, and OGLE DIA light curves 
were compared and the best one was selected in each data set as follows. First, the planetary 
deviation at IIJD-245000 = 4300-4305 was removed from each light curve, and these planet- 
free light curves were fitted with a single lens model with xallarap (binary orbital motion 
of the source). (Details are discussed in Section |3]). The photometric reduction yielding 
the smallest variance from the best model in these planet-free fits were selected to use for 
further analysis. For each data set, the error bars are rescaled so that x^/(d.o.f.) 1 in the 
planet-free single- lens fit. For CTIO V and H which have very few data points unaffected by 
the planetary deviation, this same procedure was followed including the planetary deviation 
using the best fit planetary model to all the data sets. 



Figured] shows the light curves of this event around the peak and the planetary deviation. 



3. Light Curve Modeling 

The negative deviation that triggered the initial alert is characteristic of "minor image 
perturbations", in which the image inside the Einstein ring is perturbed by a planet, and 
therefore a planet is inside the Einstein ring. In this case, two triangular caustics appear 
near the central caustic, on the opposite side of the planet, as as shown in Figure [2J The 
Danish (La Silla) data show a caustic entrance just prior to their last point, and the MOA 
and Canopus data confirm this entrance and trace its rise. From these data alone it is clear 
that the source has passed into the "depression" between the two triangular caustics and 
then passed over one of the two parallel caustic walls that bound this depressed region. See 
Figure [21 The subsequent data over the next day trace the path through a triangular caustic. 
A blind search of parameter space, in which minimizations were done with various initial 
parameters, confirms that this is the only viable topology. 
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In addition to the three single lens model parameters, the time of peak magnification to, 
Einstein radius crossing time and the minimum impact parameter uq, the standard binary 
lens model has four more parameters, the planet-host mass ratio q, projected separation d, 
the angle of the source trajectory relative to the binary lens axis a, and source radius relative 
to the Einstein radius p = ^*/6'e, or the source radius crossing time = ptE- Note that p 
can be used to estimate angular Einstein radius 9e by using t he source angular radius 9^, 
which can be estimated from its color and apparent magnitude (lYoo et al.ll2004bl ). 



A simple heuristic argument can be given to derive q, d and a from the gross characteris- 
tics of the hght curve (iGould &: Loeblll992l : iGaudi fc Gouldlll997l ). From the point-lens part 
of the light curve with the planetary perturbation excluded, we robustly find to ~ 2454311 
JD, — 55 days and Uq ~ 0.08. The time and duration of the planetary deviation is 
td — 2454303 JD and At ~ 1 day, where we adopt the duration of the negative devia- 
tion relative to the single lens model. By using these planet-model independent values, 
the position of the planet can be estimated as d_ = ( a/mj + 4 — Ud)/2 = 0.92, where 
Ud = a/ T(td)^ + = 0.166 and r(td) = {td — i^o)/^E- The angle of the source trajectory 
relative to the binary lens axis, a, can be given by sin a = Uo/ud = 0.48, therefore a = 0.5 
rad. The separation of the two triangular caustics is given by (ina.iis = 2(7 — lY^'^ in the unit 
of the planet angular Einstein radius 9p = q^^'^9-E ( Schneider et al.lll992 ). where 7 = dZ^ 
is the shear. The duration required to pass the "depression" between the two triangular 
caustics is given by At = 2(7 — l)^/^g^/^(csc a)tE- Therefore, we find that the planet has 
the sub-Saturn mass ratio, g ~ 1 x 10^^. From the light curve around JD= 2454303, we 
can roughly find the time it takes the source radius to cross the caustic tcmss — 0.23 days. 
Therefore the source radius crossing time can be estimated as t^, = tcrossSina ~ 0.1 day. 
These first order estimates of the planetary modeling are very robust. The actual light curve 
modeling will investigate several higher order effects and possible systematics, but all within 
the context of the topology defined by Figure [21 and the main conclusions remain robust. 

The light curve modeling was done by two in dependen t codes . One uses the hybrid point- 
source, individual image ray-shooting method of iBennettI ( 120091 ) . which ha s been developed 
from the first completely general finite-source binary lens calculations of iBennett fc Rhie 
( 119961 ). The other uses the same basic strategy, but was independently written by MOA. 
The bes t fit binary lens m odel was found by the Markov Chain Monte Carlo (MCMC) 
method (jVerde et al.ll2003l ). The Markov chains of preliminary runs were used to derive the 
optimal directions and step sizes for exploring parameter space. The resultant distribution 
of the chains gives us the best fit parameters and their errors. We use a linear limb-darkening 
model for the source star using the coefficients, u = 0.5250 for /-band, 0.6834 for V, 0.3434 
for H and u = 0.566 for MOA-Red which is a mean of R and /, from I Clare tl (I2OOOI ) for a 
G6 type source star with T = 5750 K and logg = 4, which is based on the best fit source 
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V — I color (see Section [5]). The best fit source and blend are plotted in the color magnitude 
diagram (CMD) (Figure |3]). The best fit standard binary lens model has a planetary mass 
ratio of g = 1.3 x 10"'* and other parameters as listed in Table [H in which q, d a and the 
source radius crossing time, t^, = pt-^ = 0.1 day are consistent with the first order estimates 
given above. 

However the overall light curve shows asymme tric residual s about the primary peak, 

which suggests either the microlensing parallax effect jGouldlll992[lAlcock et al]ll995[lsm"ith. Mao fc Wozni 
2OO2I ) by which the Earth's orbital motion distorts the light curve and/or the x allarap effect, 
which is a similar distort ion caused by the orbital motion of a binary source ( iGriest fc Hu 
I992I : iHan &: Gouldll 19971 ). Therefore, we compare the data to models that included parallax 
and xallarap. 



3.1. Microlensing Parallax 

The parallax effect is represented by two additional parameters, an amplitude tte = 
t^tci/Oe, i-e., the lens-source relative parallax vTrei = {tti — tTs) in unit of the angular Einstein 
radius 6e = Re/Di, and a direction of the relative source-lens proper motion relative to 
North toward East where Di is the distance to the lens. As shown in Table [H the best 
fit parallax model improves by Ax^ = 298 relative to the best standard binary model. If 
this parallax model were the correct model, we could derive the lens mass M = 6v,/( K7ry,) 



0.040 ± 0.005 Mq, and distance Di = AU/vr, = 867 ± 93 pc, for this model of the lens flGould 



I992I ). Here n = AG/c^AU = 8.144masMQ^ (milli-arcsec per solar mass) and we have 



assumed the source distance Dg = AU/tTs = 8.0 ± 1.4 kpc where the error is based on 17% 



standard deviation in the Galactic bar model (IHan fc Gould II2003I ). This model implies that 



the lens is a nearby brown dwarf. However, as shown in the next section, the best xallarap 
model yields a significantly better x^, with an improvement by Ax^ = 89.4. Furthermore, 
if the signal were due to parallax, we should have found the best xallarap model with the 
same (R.A.^, Dec.^) values as the celestial coordinates of the source as seen from the Earth 
(RA, Dec.)= (269.1°, -32.2°) when its period of the source orbital motion, eccentricity, and 
perihelion, celestial pole are fixed as the values of Earth's orbit. However, we obtained the 
best model with (R.A.^, Dec.^ ) = (309.4°, -24.0°) ± (0.5°, 0.2°), which is inconsistent with 
the expected values for parallax. We conclude that this distortion is not likely due solely to 
parallax. 
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3.2. Xallarap 



If the orbit is assumed circular, and the companion assumed to generate neghgible flux 
compared to the source, the xallarap effect can be represented by five additional parameters, 
an amplitude, = cis/^E, that is the semi-major axis of the source's orbit, a^, in the unit 
of the Einstein radius projected on the source plane, = ReDs/Di, the direction of the 
relative source-lens proper motion, 0^, the direction of observer relative to the source orbital 
axis, R.A.^ and Dec.^, orbital period, P^. For an elliptical orbit, two additional parameters 
are required, the orbital eccentricity, e and time of perihelion, tperi- 

The best xallarap model models, with e fixed at e = and with e as a free parameter, 
improved by Ax^ = 74 and 89, respectively, relative to the the best parallax model. 
The best fit parameters are listed in Table [H We also fitted models with a bright binary 
companion, but in every case, the dark binary companion provided the best fit. Therefore 
we keep only models having companions with negligible flux compared to the source in the 
following analysis, which would be appropriate for a white dwarf or M dwarf companion. 

In Figure m we show the of best xallarap models as a function of with orbital 
eccentricity fixed at e = and fitting for e. One can see that xallarap models are signifi- 
cantly better than the best parallax model, and that the xallarap model in which e is a free 
parameter is slightly better than the model fixing e = 0. For the xallarap models, is flat 
for > ISOdays, in which regime P^ and are strongly degenerate. 

Of course, every microlensing event must have a non-zero microlensing parallax, but the 
addition of parallax to these xallarap models did not provide a significant improvement. 
The parallax and xallarap parameters are highly degenerate and tend to complicate the 
analysis, so we have excluded parallax from most of our xallarap models. 



In Section 13. 2[ the model with the lowest is the xallarap model with non-zero e. 
However, this model leads to a xallarap amplitude of ^e = 0.35, which is larger than would 
be induced by a "normal" main-sequence companion, where ^e is expressed, making use of 
Kepler's third law, by 



3.3. Xallarap with the Kepler Constraint 




lAU / Mc 



TE \Mq 




Me 
Mc + Ms lyr 




(1) 
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From this equation and parameters for this model, the lower limit of the companion mass to 
the source is given by, 

-3^3 ( Pi. 



M, > iiri (^^J ~ SOMq, (2) 

which would imply a black hole companion, since a 50 Mq star would exceed our upper limit 
on the brightness of a companion to the source by more than 5 magnitudes. 

Black holes are quite rare compared to stars, so we should consider the prospect of a 
more normal stellar companion. We can use Kepler's third law, the projected Einstein radius, 
te, source mass, M^, a nd the source companion mass, Mr,, to constrain the magnitude of 



the xallarap vector (IBennett et al.ll2008l : iDong et al.ll2009bl ). From Section [5l we derive 



Ms = 0.9 ± 0.1 M0 and assume a white dwarf companion = 1.0 ± QAMq (which would 
be the most massive dark companion with a plausible a priori probability). Inserting these 
masses and other relevant parameters into Eq. ([T]), the maximum allowed ^e for the best 
xallarap models for the circular orbit xallarap and non-zero e models are given by 

^E,max = 0.11 ± 0.04 (e = 0) and 0.06 ±0.02 (e free), (3) 

where the error is estimated from the errors in 6'^., Mg and M^. 

Because our best fit values of ^e = 1-73 and 0.35 for the circular orbit and non-zero e 
are much larger than ^E,max given above, they are inconsistent with our upper limit on the 
source companion mass. To find the best xallarap model that is allowed by Kepler's third 
law, we have done MCMC runs with an additional constraint contribution to given by 

= ©(^E - eE.max) ' . (4) 

where ^E,max is evaluated by Eq. ([T]) with parameters in each step of the MCMC and fixed 
values of Mg = M^. = IMq and 50% error in ^E.max, which depend only weakly on other 
parameters. Here, is the Heaviside step function. 

In Table [H we show the best fit model parameters and for the circular orbit and 
non-zero e cases. In Figure H] we show the of the best-fit xallarap models with the Kepler 
constraint as a function of P^. One can see that if we impose the Kepler constraint, the 
xallarap solution with e free is better than the case of fixed e = by Ax^ = 39. Although 
this is worse, by Ax^ = 18, than the model without the Kepler constraint, this is the best 
model that is allowed for a plausible companion mass. 
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4. The Errors in Parameters with Systematics 

We have investigated second order effects to explain the clear asymmetry about the peak 
in Section I3.H 13.21 and 13.31 We also searched for models with an additional mass besides the 
lens star and planet, but neither an additional stellar or planetary companion to the lens 
star could account for observed light curve asymmetry. 

We are sure that there is an asymmetry in the light curve, because we see qualitatively 
similar trends in both MOA and OGLE light curves as shown in Figure [5l which have 
different typical seeing and were reduced by independent pipelines. However we are not fully 
confident that this xallarap amplitude is correct because of the unphysically large and an 
additional factor: there is a bright red clump giant (RCG) star with {V — I, /) = (2.01,15.56) 
at the North-East corner of the Keck AO image (see details in Section |6]) in Figure El which is 
only 1.1" away from the source. The wing of the giant star PSF interferes with photometering 
the source on the OGLE images, with typical seeing of ~ 1.2", and even worse in the MOA 
images, with typical seeing of ~ 2.0". The differential atmospheric extinction and refraction 
may cause systematic asymmetry on the light curve. Here the differential refraction causes 
the positional shift of the target on the sky, which generates residuals on the subtracted 
images in DIA. These effects depend on the color of stars and airmass. The mean airmass 
changes slowly during the event because the mean elevation of the target changes in season. 
They can generally be reduced by choosing the reference stars with the similar color as the 
target for aligning the image coordinates and solving the kernel in DIA. However, the effects 
due to the blending star with different color from the target are hard to remove. Especially, 
subtle effects on the bright blending star can cause significant effects on the faint target. So 
the photometry of this event is challenging. We tested the modeling with data points taken 
at airmass > 1.3 removed, but this does not make any significant difference. When we model 
by first removing either the OGLE or MOA dataset, the results are qualitatively unchanged. 

As argued in Section [3l the planetary deviation is clearly detected, and the planet 
parameters can be estimated robustly by simple inspection. Although our analysis of the 
parallax and xallarap fits indicate the presence of unrecognized systematic errors in the 
data, these errors do not affect these basic inferences about the planet itself. Therefore, we 
are only interested in robustly estimating the parallax or xallarap parameters to the extent 
that they can provide additional information about the lens. However, having discarded 
the parallax model for the asymmetry, the xallarap parameters themselves provide no new 
information about the lens, and are therefore of no intrinsic interest to us. We therefore 
do not further investigate the cause of the systematics in the light curve, and instead seek 
only to determine the (relatively minor) extent to which these systematics affect our precise 
determination of the planetary parameters. To do this, we consider the standard and all 
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the possible xallarap models shown in Table [T] as viable, and take the differences of the 
parameters as the size of the systematic errors. We take parameters of the xallarap model 
with non-zero e and the Kepler constraint (indicated as Xallarap^ in Tabled]), as the median. 
The resultant systematic errors are listed in Table [T] and /s,ogle = 19.51 ± 0.03 mag and 
Jb,oGLE = 18.77 ± 0.02 mag. The values and errors in the following sections are estimated 
by taking these systematic errors into account. 



5. Source Star Characterization 



We must determine the source star angular radius, in order to determine the angular 
Einstein radius, 9e, from the light curve parameters. Since we do not have infrared light 
curve data of high e nough quality to acc urately measure the source bright ness in the i nfrared 
JCould et allbood : ISennett et aPboOol l. we use the two-filter method of IYoo et all J2004bh 
to determine 9^,. The values and errors stated in this section are the final values including 
systematic errors, as they are subsequently re-estimated after Section HI However, these are 
qualitatively the same within the errors as the original results used for the initial light curve 
modehng in Section [31 



5.1. Extinction Correction 



The V and / magnitudes of the source star from the light curve fit (see Section [3]) must 
be corrected for the extinction and reddening due to the interstellar dust to infe r the spectra l 
type of the source. Because this field is out of the OGLE-II extinction map (ISumi 1120041 ). 
we estimate the extinction and redd ening to the source by using RCG, which are known to 
be an approximate standard candle (IStanek et al.lll997l : iPaczyhski fc Staneklll998l ). Figure [3] 
shows the calibrated OGLE CMD in 3x3 arcmin^ field around the event [(/, h) = (358.3488°, 
-3.6861°)]. From this CMD, we find the RCG centroid; 



(y-/,/)RC,obs = (2.14,15.70), 



(5) 



where the errors are negligible comparing to the intrinsic error in the RCG centroid, as 
described below. 

We adopt the intrinsic RCG m agnitude M^ rc n = - 0.25±0.05, My-Rn n = 0.79±0.08 and 
coloT (V — I )rc,,() = 1.04=fc0.08 f romlBennett et al.l (120081 ) which is based on lCirardi fc Salaris 



(I2OOII ) and ISalaris fc Girardi I fl2002h . where the error is assigned based on the size of the 



theoretical corrections to the RCG magnitudes. Taking account of the bar structure of the 
Galactic bulge, the offset of the distance modulus (DM) between the GC that is assumed to 
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be at 8.0 ± 0.5kDC (lReidl ll993r) and the average stars in the field have ADM = 0.00 ± 0.05 
( INishiyama et aLll2005[ ). So the dereddened RCG centroid in the field is expected to be 



{V - /, /)rc,o = (1.04, 14.27) ± (0.08, 0.15). 



(6) 



Comparing these centroids (Eqs. [5] and [6]) , we find the average reddening and extinction 
in this field is 

{E{V - I),Aj) = (1.10, 1.43) ± (0.08, 0.15), (7) 

wher e Ryi = Ay/ E(y —I) = 2.30, which corresponds to i^y = 2.64 (ICardelli. Clayton &: Mathis 
19891 ). Applying this average extinction to this event, the source's {V — I, /)s,ogle from fit- 
ting of the well calibrated OGLE V and / light curve and the dereddened source magnitude 
and color {V — I, /)s,o are 



(V^-/,/)s,OGLE 
s,OGLE,0 



(1.85, 19.51) ±(0.06, 0.03), 
(0.75,18.08) ± (0.10,0.16), 



(9) 



Independently, the dereddened source color, (V^ — /)s,ctio,o = 0.77 ±0.02 ±0.08, is estimated 
by comparing (V — /)rc,o; the CTIO RCG color and the CTIO source color {V — /)s,ctio 
which is given by the model-independent regression of CTIO V and / light curves. This is 
consistent with {V — /)s,ogle,05 but more accurate. In the following analysis, we adopt the 
value: 

{V - I, /)s,o = (0.77, 18.08) ± (0.08, 0.16). (10) 



which implies that th e source is a mid-G st ar in the bulge ( iBessell &: Brettlll988l ) with mass 
of M, = O.9±O.1M0 f |Schmidt-Kal^ll982[ ). The reddened {V-I,I)s is plotted in Figure El 



The dereddened blended light from the light curve is 

(V, /)b,o = (17.71, 17.34) ± (0.18, 0.15). 



:ii) 



Note that if this blended light is from the lens or companion of the lens, these values may 
be over-corrected for extinction because these objects are in front of the source. Thus, these 
magnitudes can be used as an upper limit on the combined light of the lens, any companion 
of the lens, and the source in the following analysis. 



5.2. Source Star Angular Radius 



Following IYoo et al.l ( l2004bl ). the dereddened source color and brightness {V-K, -ft')s,o = 
(1.69,17.16) are estimated using the observed {V-I, J)s,o as given by the equation ffTOj) 
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and the color-color relation (IBessell &: Brett 



1988h. By u s ing th is {V-K, K\^o and the 



empirical color/brightness-radius relation by iKervella et al.l (|2004[ ). we estimate the source 
angular radius, 6^ = 0.81 ± 0.07 fias, where the error includes uncertainty in the color 
conversion and the color /brightness-radius relation. On the other h and, {V-I, /)s,o and 
optical color/brightness-radius relation by IKervella fc Fouqua (|2008[ ) yields 6** = 0.83 ± 



0.05 /ios, which is consistent with above. We adopt the mean of these estimates. 



= 0.82 ±0.07/ias. 



(12) 



The angular and projected Einstein radii, and lens-source relative proper motion fi are 
estimated, respectively, as: 



Oe = = 529 ±84/ias, (13) 

TE = 9exDs = [4.2 ± 0.7] (^^^ AU. (14) 

/i = — = 3.5 ± 0.6 masyr" . (15) 

As shown in the top right panel of Figure U\ the measured value of /i = 3.5 masyr^^ is 
typical for the bulge lenses but smaller than the typical value for disk lenses, 5-10 masyr"^, 
although it is not inconsistent with a disk lens. 



6. Keck AO Observation 

H and K AO images of the event were taken by the Keck telescope at HJD=2454332. 77689 
and 2454332.77977, respectively. The magnification at the time of the Keck images are taken 
is ^Keck = 2.490. The magnified source position on the OGLE difference image is marked 
with the error of ~2.5 pixels (25 mas) in i^'-band. From the Keck K-band image, the density 
of ambient stars with 3-0" detection limit that correspond to K <18.1 mag, is ~0.3 arcsec"^. 
We conservatively assume that the separation of two stars must be more than the measured 
FWHM of the PSF of 0.08 arcsec in order to be separately resolved. Therefore the prob- 
ability of blending with any random interloper, that is not related with this event, is only 
~0.6%, implying this object is almost certainly the source, the lens and/or their compan- 
ion. The H and K magnitude were measured by PSF photometry and calibrated to 2MASS 
system using t he H and K i r aages taken by the IRSF telescope in South Africa, following 



the method in iJanczak et al.l (120101 ) 
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{H, i^)s,Keck = (16.53, 16.23) ± (0.03, 0.02), (16) 

and the magnitudes corrected for extinction given by Ah I Ay = 0.176 and AkI Ay = 0.105, 
which are estimated by using ICardelh. Clayton fc MathisI (1l989[ l's law with the Ry = 2.64 
measured above, are 



{H, ir),,o,Keck = (16.09, 15.96) ± (0.04, 0.03). 



(17) 



The I — H and I — K source colors are estima ted from {V — / ) ., n giv en by light curve 
fitting (Eq. [10]) by using the color-color relations of iBessell &: BrettI ( 119881 ) . 



{i-H,i~ = met'oil o.92l°:}2). 

Therefore H and K source magnitudes are given as, 

{H, = (17.231°;?°, 17.16 ± 0.20). 
Then, the magnitude of the source when the Keck images were taken are 

(H, ir),,o - 2.5(log[AKeck], log[AKeck]) = (IQ-^St'oH 16-17 ± 0.20) 



:i^ 



(19) 



(20) 



By subtracting Eq. f l2U]) from Eq. f lT7|) . we have the magnitude of the lens and/or 
companion of the lens or source, which serve as an upper mass limit of the lens 



max,0 



18.31^9, 17.9lo~ ) 



(21) 



This i^-band magnitude implies t hat the upper limit of the lens is ari early G dwarf with 
mass of M,,max = I.OI^^Mq from JSchmidt-Kalei]ll982l : lBessell fc Brettlll988h . These H and 
i^-band upper limits are used for constraining lens star in Section [71 If we could obtain a 
second epoch AO observation that gave us the baseline photometry, we would be able to 
constrain {H, -ft')i,max,05 much better. 

For other (brighter) events, we have found that the if-band source magnitude estimated 
by fitting the CTIO i7-band light curve gives a more precise value for the H magnitude of the 
source. But, when we attempt such an analysis for this event, we find significant indications 
of systematic errors. This is not very surprising because this target does not reach high 
magnification and is heavily blended with a nearby bright star. Also, because of the bright 
infrared sky brightness, the CTIO if-band images do not go as deep as the optical images. 
Therefore we do not use this CTIO if-band source magnitude in the following analysis. 
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7. Lens System Masses and Distance 



The lens system mass, M;, distance, D/, and lens-source relative velocity are directly 
constrained by only two measured parameters, the Einstein radius crossing time, ^e, and the 
angular Einstein radius, 9e- However we can further constrain them by a Bayesian analysis 
using a model of Galactic kinematics (lAlcock et al.lll995t iBeaulieu et al.ll2006l : iGould et al. 
20061 : iBennett et al.ll2008l ). The mass of the planet can be determined to the same precision 
as Ml because the uncertainty in the mass ratio, g, is much smaller than the uncertainty in 
the Bayesian estimate of Mi. 



7.1. Planetary System Parameter for OGLE-2007-BLG-368Lb 

For this event, we observed finite source effects from which we measured the angular 
Einstein radius 6'e (Eq. [13]), or equivalently the proper motion /i (Eq. [15]), of the lens system. 
So we can break one link of the three-fold degeneracy by the relation, 

e% = KMvTrel. (22) 

To produce the likelihood distributions shown in Figure U\ we compute the likelihood 
by combining this eq uation and the measured values of 6e and with the Galactic model 
(IHan fc Gould I l2003l ) assuming the distance to the GC is 8 kpc. Here systematic errors 
in parameters estimated in Section H] are taken into account. This analysis yields that the 
primary is a K-dwarf with mass of Ml = OMtofe Mq at A = 5.9+?;^ kpc and a planetary 
mass of Mp = x 10"^ Mq = 20^8 and projected separation of r± = 2.8toi AU. 

The physical three dimensional separation a = S.S^Qg AU, car i be estimated statis tically 



by putting a planetary orbit at random inclination and phase (IGould fc Loeblll992l ). The 
lens-source relative proper motion fi = S.S^Qg masyr"^, which is consistent with the value 
given by Eq. f[T5l) , favors that the lens is in the bulge rather than the disk in which typically 
/i = 5-10 masyr"^. 



7.2. Comparison to Other Known Exoplanets 

Figures M and M compare the masses and semi-major axes of the planets found by 
microlensing to those found by other methods. Figure M takes into account the different 
masses of the primary stars and uses the ratio of the semi-major axis to the position of the 
snow-line as the x-axis parameter in order to display the data in a way more relevant to 
planet formation theories. 
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The positions of all the microlensing planets on these plots are determined by a Bayesian 
analysis similar to the one we present for OGLE-2007-BLG-368Lb. However, there is a crucial 
distinction. The events plotted with red filled circle and error bars have masses determined 
either by microlensing parallax measurements or by direct identification of the lens star in 
HST images, so these can be considered to be actual measurements. The other microlensing 
planets, plotted with red open circle and error bars, are like OGLE-2007-BLG-368Lb, in 
that the light curve measurements do not directly determine the lens system mass. For 
these events, the derived parameters have a significant dependence on the assumed prior, 
and we must be careful not to over-interpret the results. For example, we cannot use the 
results of these Bayesian analyses to study the probability that a star will host a planet as 
a function of its mass, because these estimates of the host star mass depend upon our prior 
assumptions about this planet hosting probability. Instead, such questions must be studied 
with a new Bayesian analysis using only directly measured quantities as constraints. 

There are planetary microlensing events that warrant some additional discussion. The 
Bayesian analyses for these events yield double-peaked likelihood functions. This gives results 
that are extremely sensitive to the prior assumptions, so one should not directly use the 
Bayesian results in these cases. Event MOA-2007-BLG-400 has a severe o? O l/c? model 
degeneracy, which yields a factor of 10 uncertainty in the projected star-planet separation. 
We extend the error bars from the 1-a lower limit on the semi- major axis from the d < 1 
solution to the 1-a upper limit from the d > 1 solution. 



The other ambiguous planet is MOA-2008-BLG-310Lb fiJanczak et al.ll2010r ). This event 
is unusual because the kinematics favors a low-mass primary of ~ O.IM©, while the excess 
flux seen in VLT/NACO images of the source star suggests a much more massive planetary 
host star with M ~ 0.7 Mq. But this excess flux could be due to a companion to the lens, 
source, or the chance superposition of an unrelated star. So, the Bayesian analysis yields two 
peaks for the lens star (and planet) masses, but the relative weighting of these two peaks is 
quite sensitive to the assumed prior. So, as with MOA-2007-BLG-400, we use the 1-a upper 
and lower limits on the high-mass and low- mass solutions for our error bars for this event. 
For the central point, we use the geometric mean of the peaks of the high- mass and low- mass 
solutions. 



8. Constraints on the Planetary Mass Function 

In Figures [8] and |9l it appears that the distribution of planets found by microlensing 
is roughly independent of mass above IM0, with perhaps a peak at M ~ lOM®. However, 
the probability that a planet can be detected by microlensing depends on its mass, and 
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these figures have not been corrected fo r the planet detection efficiency (lAlbrow et al.l 12000 
Gaudi &: Sackettll2000l : iRhie et al.ll2000l ). A full calculation of the planet detection efficiency 
(Cassan et al., in preparation) including detail assessments of various potential systematics 
is beyond the scope of this paper, but we can obtain interesting constraints on the planetary 
mass function using a simple model for the relative planet detection efficiency. 

For events with signals due to the planetary caustic ( Gould fc LoeblE992 ). there are some 
simple theoretical arguments regarding the dependence of the planet detection efficiency on 
the mass ratio, q. If we igno re finite source effects , which are usually unimportant for 
planets with masses ^ lOM® ( iBennett fc Rhielll996l ). then the planetary caustic shape is 
nearly independent of q, and its area scales as q. We can define a planet detection region 
as the area of th e lens magnification pattern that differs from the single lens light curve 
( iPaczyhskil Il986[ ) by more than some threshold (either relative or absolute). With such a 
definition, the area of the planet detection region will scale as the size of the planetary 
caustic, as q. Then the probability that a given source trajectory will cross the planet 
detection region scales as the linear dimension of this region, which goes as g^/^. So, in the 
limit of very good light curve coverage, the planet detection efficiency for planetary caustic 
events should scale as g^/^. However, the duration of the planetary deviation also scales 
as g^/^, and with sparse light curve coverage or large photometric error bars, the detection 
efficiency can scale as steeply as q. In general, we expect that situation to be in between these 
extremes for planetary caustic events, so that we should expect that the planet detection 
efficiency should have some scaling intermediate between g^/^ and g. We have ca l culate d 
the detection efficiency for OGLE-2007-BLG-368 using the method of iRhie et all J2000h . 
Of course, it would be inconsistent to use the follow-up data that was taken because the 
planetary signal was recognized in such a calculation, so we have only included the regularly 
scheduled survey data in this calculation. The dependence of the detection efficiency on 
the detection threshold in different alert systems is negligible compared to the dependence 
on the lightcurve coverage of the dataset. This calculation gives a detection efficiency for 
OGLE-2007-BLG-368 that scales as ~ g° ® at the range of g appropriate for Neptune-mass 
to Jupiter-mass planets. This same scaling also holds true for the two other microlensing 
planets discover ed through planeta ry caustic deviations, OGLE-2003-BLG-235, and OGLE- 
2005-BLG-390 (iKubas et al.ll2008l ). For all the calculations in this section, we assume that 
the distribution of planets is uniform in log{d) for all separations, d, and we sum over all 
separations. 

The situation is somewhat different for high magnification microlensing events, which 
are detected throu gh perturbatio ns of the central caustic. Since the linear size of the central 
caustic scales as g (lDominiklll999l ). one might expect that the detection efficiency would scale 
more steeply than g for the same reasons that the planetary caustic planet detection efficiency 
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scales more steeply than g^/^, but this is not the case. The reason for this is that for events of 
sufficiently high magnification, A^ax > 50 or so, the planet detection efficiency for Jupiter- 
mass planets saturates at 1 for a wide range of separations. This is, in fact, the main r eason 



why the observing groups focus on high magnification events (iGriest &: Safizadehlll998l ). The 



planet detection ef 


iciency has been calculated for a number of high magnification events 


(Albrow et al. 200] 


Gaudi et al. 2002: 


Rattenburv et al. 2002: Abe et al. 


2004: Yoo et al. 


2004b: Done et al. 


2006: Bennett et al. 


20081: Naeava 20091: IVee et all 200 


9: Batista et al. 


2009). and these events reveal detection efficiencv scalines that ranee from for MOA- 



2006-BLG-130 and MOA-2007-BLG-192, to q''-" for OGLE-2008- BLG-279, which is the event 



with the highest planet detection sensitivity ( lYee et al.l l2009l ). Generally, the scaling is 
shallowest for the events with the highest sensitivity to planets and steeper for events with 
lower sensitivity due to lower peak magnification, less complete light curve sampling, or less 
precise photometry. For the collection of high magnification events observed, we estimate 
that the mean detection efficiency scales as gO-5±o i^ qj^({ for all microlensing events, we 
estimate that the detection efficiency scales as g0-6±0-i. 

We can now use the detection efficiency estimate to infer some properties of the dis- 
tribution of planets in our Galaxy. In analogy to the stellar mass function, we define the 
planetary mass ratio function, dN^\/ d\ogq, such that the number of planets per star in a 
logarithmic mass ratio interval is given by dN^\/ dXogq. We assume to have a power-law form 
for the planetary mass ratio function, 

dN^x 



N^q^Q{q-qo)Q{qi-q) 



(23) 



dlogq 

where go and qi are the lower and upper limits on the planetary mass ratio, (go could alter- 
natively be considered to be a low- mass-ratio cutoff on the planetary detection efficiency.) 
Thus, n = would imply that there are an equal number of planets in every logarithmic 
mass interval, and n = —1 would imply that total mass of planets in every logarithmic mass 
interval is the same. 

We can estimate the parameters, A'^o and n that describe the planetary mass ratio func- 
tion using a likelihood analysis. The expression for the likelihood function for the planetary 
mass ratio function parameters is just the Poisson probability of finding the observed num- 
ber of events, A'obs, times the product of the probability of finding events which each of the 
observed mass ratios, g,. This can be written as 



TT 



n 



dNr, 



dlogq 



(24) 



0.6±0.1 



is the planet dete ction efficiency and Ap^p is the number of events 



where £{q) oc g' 

expected for the given Ao and n values (lAlcock et al.lll996l . 119971 ). However, since we have 
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only calculated relative and not absolute efficiencies, we cannot calculate N^xp and we cannot 
hope to constrain Nq. Therefore, we adjust $0 so that N^xp = N^hs, and evaluate the 
likelihood function for only the power-law index, n, of the planetary mass ratio function. 
The resulting likelihood function based on the ten planets discovered by microlensing is shown 
in Figure [TOl and the resulting planetary mass ratio function index is n = —0.68 ±0.20, with 
a 95% confidence level upper limit of n < —0.35. The core of this distribution is similar to 
a Gaussian, but the distribution is skewed, with a higher probability of a > 2-a deviation 
at small n than at large n This error bar includes the ±0.1 uncertainty in the detection 
efficiency power law index {£{q) oc gO-6±o i) This result does have some dependence on our 
choice of the lower and upper cutoffs of go = 3 x 10^^ and qi = 0.015, but the variation due 
to the choice of these cutoffs is much smaller than the resulting uncertainty in n. 

This result for the power law index indicates that we should expect times as many 
cold Neptunes (g ~ 5 x 10~^) as Jupiters (g ~ 10"'^), with a 95% confidence level lower limit 
of 2.8 times as many cold Neptunes as Jupiters. This is in line with the basic pred ictions of 
the core accretion model (jida fc Linll2004l : iLaughlin. Bodenheimer fc Adamsll2004l ). as these 
models predict a large population of Neptune-like, "failed Jupiter" cores to form beyond 
the snow line, particularly for stars of less than a solar mass, which make up most of the 
sample probed by microlensing. However, it still may be po ssible to explain this result in 
the context of the gravitational instability theory (jBossll2006l ). 



This power law index of n = —0.68 ±0.20 is ste eper than (but consistent with) the index 
of n = — 0.31±0.20 found by lCumming et al.l (120081 ) for more massive planet s orbiting mostly 



solar-t ype stars. This is also steeper than the mass function prediction of iMordasini et al. 



(I2OO9I ) for solar-type stars, although this theoretical mass function is not a power law. Radial 
velocity surveys also find that hot Neptunes, with periods less than 50 days are quite common 
around G and K dwarfs (IMayor et al.ll2009l ). 



9. Discussion and Conclusions 

We have presented the analysis of the OGLE-2007-BLG-368 planetary microlensing 
event, which indicates that the planet OGLE-2007-BLG-368Lb is a Neptune-mass planet. 
We also find evidence for low level systematic errors in the light curve, which however do not 
affect this conclusion. We estimate the systematic errors by taking the differences between 
the various models, i.e., the standard and xallarap with and without Kepler constraint. 

By using a Bayesian analysis, we found the planet has a mass of Mp = 20^g and 
a projected separation of r± = 2.8lo;g around a K-dwarf with mass of Mi = 0.641q;26 
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Mq at Di = 5.9tiA kpc. This is the 4th Neptune -mass planet detected by microlensing. 
In Figure El we plot these planets as a function of mass vs. semi-major axis along with all 
known exoplanets. FigurelHlis the same as Figure|H]but the semi-major axis is divided by the 
snow line, which is taken to be at Osnow = 2.7 AU M/Mq. As for the microlensing planets in 
this figure, we are starting to see a broad concentration of ~ 10 planets beyond the snow 
line. This is as expected from the core accretion theory. This theory predicts that the most 
massive solid planetary cores should form beyond the snow line, which then accrete nebular 
gas and become the gas giants around the solar-type star. On the other hand, they become 
Earth-mass to Neptune-mass icy rocky planets around M-dwarfs. Comparing four Neptune- 
mass, five Jovian planets and one between Neptune and Saturn found by microlensing, it 
confirms that col d Nep tunes are relatively common around low mass primary stars analyzed 



bv lGould et al.l (120061 ). 



We have presented an analysis of the exoplanet mass ratio function. Assuming that the 
number of planets scales as a power law in the mass ratio, q, we define the mass ratio function 
as dNpi/dlogq oc g" power law mass function, over the mass range of a few Jupiter masses 
down to a few Earth-masses, we find a power law index of n = —0.68 ± 0.20, which indicates 
that Neptune-mass planets are substantially more common than Jupiter-mass planets. 

The planetary signature of this event was detected in real-time in the data points from 
survey telescopes MOA and OGLE. Then the signature was greatly clarified by intensive 



follow-up observations prompted by the alert. T 
the second of its kind after OGLE 2005-BLG-390 ( 



lis is a planetary ca ustic crossing event. 



Beaulieu et al.ll2006[) among all planetary 



microlensing events. OGLE 2003-BLG-235 (IBond et al.l |2004| ) crossed the planetary part 
of a resonant caustic. Although the time of the planetary deviation in these events can 
not be predicted for planetary-caustic events, the potential event ra te is higher than c entral- 
caustic event in which the time of the planetary deviation is known (IHan fc Kimll200ll ). This 
discovery shows that the high cadence survey observations that MOA is conducting, have 
a great potential to increase the event rate of the planetary microlensing. In 2010, OGLE 
will upgrade its camera to 1.4 deg^ FOV (OGLE-IV), which will enable OGLE to follow a 
similar strategy of high cadence monitoring for planetary signals. 

Multi-continent high-cadence observing will commence in 2010 with the start of the 
OGLE-IV project, and in future years will expand further when the Korean Microlensing 
Telescope Network (KMTNet) is commissioned. These improvements can be expected to 
dramatically increase the number of microlensing planets, and in particular those like OGLE- 
2007-BLG-368Lb, that are discovered via planetary-caustic perturbations. 
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Fig. 1. — The light curve of OGLE-2007-BLG-368 over the whole event (top-panel), around 
the planetary deviation (lower-left panel) and the second caustic crossing (lower-right panel) 
with the residual from the best fit model. The red lines indicate the best fit xallarap model 
with the Kepler constraint (see Section [3^ . Here the light curves of /zFUN CTIO /, H and 
PLANET Brazil are binned by 0.01, 0.01 and 0.02 days, respectively, for clarity. Note that 
the fittings were carried for un-binned light curves. 
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Fig. 2. — The caustics (red lines) and critical curves (black lines) of OGLE-2007-BLG-368 
for the best model fitting e with the Kepler constraint (see Section I3.3p . The blue line 
represents the trajectory of the source. The inset shows the zoom around the planetary 
caustic crossing, where the gray scale indicates the magnification pattern. The circle in the 
inset represents the best fit source size. 
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Fig. 3.— The OGLE {V-I, I) color magnitude diagram around OGLE-2007-BLG-368. The 
filled circle and triangle indicate the source and blended light from the fit, respectively. The 
filled square represent the total fiux of the source and blend. Here the errors in / are too 
small to be visible. The cross indicates the center of the RCG. 
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Fig. 4. — The of the best xallarap model as a function of the orbital period of the source 
star and its companion. The solid and dashed lines indicate the model with fixed orbital 

eccentricity e = with and without the Kepler constraint, respectively. The dot-dashed and 
dotted lines indicate the model allowing a free-fit of e subject and not subject to the Kepler 
constraint, respectively. The best parallax model is plotted as a for comparison. 
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Fig. 5. — The light curves of MOA Red (top panel) and OGLE I (bottom panel) with the 
best standard (red line and residual) and xallarap models with e being fit subject to the 
Kepler constraint (blue line and residual). Here MOA data are binned by 1 day outside of 
the planetary signal at JD-2450000=4300-4304. We can see a similar asymmetry in both 
light curves is well fitted by the xallarap model in both cases. 
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^ 10.2 X 10.2 " narrow cam 

Fig. 6. — K-band Keck AO narrow camera image of OGLE-2007-BLG-368. The magnified 
source position on the OGLE difference image is marked as "T" with error of ~2.5 pixels 
(25mas) in i^-band, where the coordinate is ahgned by using the 5 brightest stars in the 
K-band image. This object is almost certainly the source, the lens and/or their companion. 
Blending with random interloper is unlikely with this stellar density (see Section [6]). The 
bright RCG star marked as "1" is 1.1" way from the source, whose PSF tail covers the source 
on OGLE image with typical seeing of 1.2". 
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Fig. 7. — Probability distributions from a Baycsian analysis for the distance, Di, transverse 
velocity, Vt, mass, Mgtar, Einstein radius (-Re), V, I, H and iT-band magnitudes of the 
primary star of the lens system. The vertical solid lines indicate the median values. The 
dark and light gray shaded regions indicate the 1-cr and 2-a limits. The gray solid and 
dashed curves in the top-left panel indicate the mass-distance relation of the lens from the 
measurement of 9-e with l-cr errors, respectively, assuming Dg = S kpc. Note Dg is not fixed 
in the actual Bayesian analysis. Thick solid and dashed lines in the top-right panel represent 
the typical /i distributions of the bulge and disk lens populations, respectively. The vertical 
dashed and dotted lines in the V, I, H, and K-hand panels represent observed upper limit 
and 1-0" error, respectively. 
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Fig. 8. — Known exoplanets as a function of mass vs. semi-major axis, along with the 
predicted sensitivity curves for various methods. The red filled and open circles with er- 
ror bars indicate the microlensing planets with mass measurements and mass estimated by 
Bayesian analysis, respectively (see section \T72^ . OGLE-2007-BLG-368Lb is indicated by the 
gold-filled open circle. The blue dots represent the planets first detected via transit. The 
black bars with upward-pointing error bars (indicating 1 a sin i uncertainty) are the radial 
velocity planet detections. The green and magenta triangles indicate the planets found by 
timing (including the pulsar planets) and by direct detection, respectively. The yellow, cyan, 
and light green shaded regions indicate the expected sensitivity limits of the radial velocity, 
Kepler and SIM space missions. The red and pin k curves indicate the p redicted lower sensi- 
tivity limits for a ground-based and space-based (IBennett fc Rhidl2002l ) microlensing planet 
search program, respectively. The solar system's planets are indicated with black letters. 



-37- 



m 
03 



Exoplanet Discoveries vs. Snow-Line 
iiii^ ■ 



103 ^ 



CO 
CD 

m 102 

s 

u 
cd 



10 E 



1 - 



0.1 E 



1 



1 1 I I I I I I li -L , I |J, ll-J-lU llljt I 1 11 I I 1 I I ' I I M 



1 ill 



JJUJ 




J 



J I 



0.01 0.1 1 10 

(semi-major axis)/(snow-line) 

Fig. 9. — Known exoplanets as a function of mass vs. semi-major axis divided by the snow 
line, which is taken to be at Osnow = 2.7 AU M/Mq. As in FigureEl microlensing planets are 
indicated by red filled and open circles with error bars (see section . OGLE-2007-BLG- 
368Lb is indicated by the gold-filled open circle. Blue dots represent the planets first detected 
by transits. The black bars with upward-pointing error bars are the planets detected via 
the radial velocity. The green and magenta triangles indicate the planets found by timing 
(including the pulsar planets) and by direct detection, respectively. 
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Fig. 10. — The probability distribution of the power law index, n, of the planetary mass ratio 
function, "^{q), based upon the mass ratios of the ten exoplanets detected by micro lensing and 
our estimate of the planetary detection efficiency. This calculation yields n = —0.68 ± 0.20, 
with a 95% confidence level upper limit of n < —0.35. 
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Note. — Standard model includes neither parallax nor xallarap effects. HJD' = HJD-2450000. Models with super- 
scripts "*" and "ii"" indicate fixed e = and Kepler constraint Mg = = IMq, respectively. The lines with "a" 
list the 1-0" error of parameters given by MCMC, for which the xallarap parameters are fixed at the best values for the 
xallarap models because xallarap parameters are strongly degenerate and it is hard to satisfy the convergence criteria, 
^systematic iudicatcs the systematic errors (see Section H]). 



